Distribution of the Riemann Zeros Represented by the Fermi Gas 
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The multiparticle density matrices for degenerate, ideal Fermi gas system in any dimen- 
sion are calculated. The results are expressed as a determinant form, in which a correlation 
kernel plays a vital role. Interestingly, the correlation structure of one-dimensional Fermi gas 
system is essentially equivalent to that observed for the eigenvalue distribution of random 
unitary matrices, and thus to that conjectured for the distribution of the non-trivial zeros 
of the Riemann zeta function. Implications of the present findings are discussed briefly. 
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In this study we consider the noninteracting A-particle fermion system at zero 
temperature in any spatial dimension d. In the following analysis the spin-polarized 
case is considered, and the volume of the system is taken to be unity for notational 
simplicity. The A-body density matrix-^ for this Fermi system is then expressed in 
terms of the product of the Slater determinants with the plane waves^'^ as 
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where rj and ki refer to the spatial coordinates and the wavenumber vectors, respec- 
tively; 
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denotes their signature. The density matrix given by Eq. (1) satisfies the normal- 
ization condition as 
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Here, noting the identities for N x N matrices A, its transposition A T , and -B 
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is an N x iV matrix whose components are 
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for 1 < i, j < N . The n-body (1 < n < iV) density matrix^* is then given by 
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with an n x n matrix K^ n \ It is noted that the prefactor before det.K^ n ) in Eq. 
(10) becomes unity for finite n and ./V — > oo. 

A problem in the following is to find an explicit expression for the correlation 
kernel, 
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in any dimension d, which represents Eq. (9) in the limit of infinite N. Here, the k 
summation is expressed by 
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and for the degenerate Fermi gas, 



n k = 9(k F -\k\) (13) 

refers to the occupation number density in the wavenumber space with 9{x) and k F 
being the step function and the Fermi wavenumber, respectively. 

The occupation number density satisfies the normalization condition, 
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Therefore, the Fermi wavenumber k F is given by 
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is the volume of the d-dimensional sphere with the radius of k F ; r(s) is the gamma 
function.®'® The Fermi wavenumber is thus expressed as 
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By choosing the direction of r in parallel with the d-th component of the 
wavenumber vector k, Eq. (11) is expressed as 
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we find an explicit expression for the correlation kernel as 

K(r;v) = r(v + l)(JtS\ J u (k F r) (21) 

with v = d/2, where J v {z) is the Bessel function of the first kind of order i/P® 
It is remarked that, though Eq. (21) has been derived for integral values of spatial 
dimension d, K(r; v) may be regarded as a continuous function of the auxiliary order 
variable v, which follows from its series representation.®'® 

Let us here consider the case of d = 3 and v = 3/2. The correlation kernel then 
reads 
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is the spherical Bessel function of the first kind of order nP® This expression, Eq. 
(22), combined with Eq. (10) in the case of two-body density matrix, reproduces 
a well-known resuhP''^ 1 for the pair distribution function of the three-dimensional 
Fermi gas. 

In the case of d = 1 and v = 1/2, we find 
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Recalling Eq. (10), this correlation structure is essentially the same as that for the 
eigenvalues of random matrices in the Gaussian unitary ensembleP 1- ® Interestingly, 
it has been known that this type of correlation structure may hold also for the 
distribution of zeros in the Riemann zeta functionP 

The Riemann zeta function^ for complex variable s is defined by 
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for Re s > 1, where n and p mean the natural numbers and the prime numbers, 
respectively. After the analytic continuation over the whole complex plane, the 
has non-trivial zeros in the critical strip, < Re s < 1, and the Riemann hypothesis 81 
states that all of them lie on the critical line Re s = 1/2; that is, 

cQ+i*)=0 (26) 

has non-trivial solutions only when t = tj (j = 1, 2, • • • ) are real. 

The mean density of the non-trivial zeros of £(s) increases logarithmically with 
height t up to the critical line. We then define unfolded zeros by 

Wi = |log|. (27) 

It has then been conjecturecP'^ 1 that the pair correlation function of the unfolded 
zeros for j — >■ oo has a form, 

2 



„(o\ , n / sin TTW \ 

rM(w) = 1-( , 28 

TTW 



which is analogous to that for the one-dimensional (d = 1) Fermi gas, that is, Eqs. 
(24) and (10) for n = 2. This conjecture has also been generalized for all the n-point 
correlations as a determinant form analogous to Eq. (10) P^'^ 

Thus, the case of v = d/2 — > 1/2 in the Fermi gas system gives a special 
correlation structure observed in the random unitary matrices and the Riemann zeta 
function. Therefore, the behaviors of the multiparticle correlations of the Fermi gas 
system with the correlation kernel K (r; v) at and around v = 1/2 may provide useful 
information about the random matrices and the zeta function, regarding K (r; v) as 
a continuous function of v. Another challenge, which would be more ambitious, 
is to look for a family of functions whose zero distributions are described by the 
correlation functions given by Eqs. (10) and (21) for arbitrary values of v. 

In summary, it has been revealed in this study that the correlation structure 
similar to those observed in the eigenvalues of the random unitary matrices and in 
the Riemann zeros is embedded in the Fermi gas system as well. 
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